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Abstract 


The  work  reported  here,  for  the  period  12/5/84  through  4/30/86,  produced 
the  following  results: 

(a)  A  theoretical  formulation  of  the  flow  of  plasma  in  a  variable  area 
accelerator  under  conditions  where  the  voltage  is  dominated  by  the  back 
e.n.f. ,  showing  novel  features  akin  to  those  found  in  ordinary  gas 
dynamics,  but  with  the  aagnetoacoustic  speed  playing  the  controlling  role. 

(b)  A  numerical  model  of  an  axisysmetric  MPD  thruster  of  realistic 
geometry  with  fully  coupled  gas  and  electrodynamic  effects,  but  limited  by 
numerical  difficulties  to  conditions  well  below  onset. 

(c)  Design  and  construction  of  test  channels  to  investigate  the  effects 
predicted  by  the  above  theories,  and 

(d)  Generation  of  a  limited  computerized  MPD  data  base. 
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l.  Introduction.  Research  Ob  iectives 

Our  work  statement  for  the  third  year  of  our  AFOSR  Grant, 
covering  the  period  of  Dec.  15,  1984  to  Dec.  14,  1985  (later 
extended  to  April  30,  1986)  included  three  broadly  defined  tasks: 

a)  .Analytical  study  of  simplified  models  of  the  flow  and 
current  attachment  in  MPD  channels,  with  emphasis  on  the  low  density 
regions  of  the  lip  anode  3nd  the  plume. 

(b)  Development  of  a  2-D  axisymmetric  numerical  model  for 
calculation  of  the  steady  state  flow  and  electrodynamics  of  realistic 
MPD  thruster  configurations. 
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(c)  Experimental  testing  of  theoretical  predictions  on  MIT  designed 
thrusters  to  be  operated  at  the  E&D  Associates  facility  in  Alexandria, 
Virginia. 

During  the  year,  a  fourth  task  was  added,  to  construct  a  Personal 
Computer  based  Data  Base  of  MPD  data  and  design  information. 

In  the  following  sections  we  give  an  overview  of  the 
accomplishments  in  these  tasks,  making  reference  to  several  Appendices, 
where  additional  details  are  contained. 

2-  Outline  of  Work  Accomplished 

Early  in  our  efforts  to  develop  an  understanding  of  the  way 
current  distributes  itself  in  an  MPD  thruster,  and  the  way  mass  tends 
to  migrate  towards  the  central  cathode  and  leave  partially  depleted 
anodic  regions  behind,  it  became  apparent  that  the  role  of  magnetic 
convection  as  a  dominant  effect  in  the  MPD  regime,  while  implicitly 
acknowledged  in  previous  work,  had  not  been  properly  appreciated, 
nor  had  the  simplifications  afforded  by  this  dominance  been  adequately 
explored  before.  In  essence,  this  indicated  the  desirability  of 
examining  the  high  Magnetic  Reynolds  Number  regime,  which,  together 
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with  the  dominance  of  aagnetic  over  fas  preaaura,  can  be  seid  to 
defiae  the  specific  MPD  reals.  The  practical  value  of  this  Magnetic 
Reynolds  Rusher,  Ra,  is  of  the  order  of  3  to  10,  which  is  only 
soderately  high,  but  the  clarifying  insights  afforded  by  the  high 
Ra  liait  are  still  expected  to  be  quite  useful  in  steering  future 
research  and  design  work. 

A  first  result  of  this  redirection  of  our  analytical  work 
was  the  AIAA  paper  included  here  as  Appendix  1.  In  it  we  first 
show  how  the  Rn->  •  liait  reduces  the  plasaa  dynasics  to  a  fora 
fasiliar  froa  ideal  gas  dynasics,  since  aagnetic  pressure  is 
found  to  be  governed  by  density  alone.  We  thus  are  led  naturally 
to  concepts  such  as  choking  (but  at  the  Alfven  velocity  instead 
of  the  speed  of  sound),  convergent-divergent  accelerators,  with 
the  lofting  being  dictated  by  the  desired  current  distribution, 
and  anode  lip  current  concentrations  as  a  fora  of  Prandtl-Meyer 
expansion  fans.  Further,  to  explore  the  two-diaenional  effects 
induced  by  Rail  currents,  we  then  introduce  perturbations  of 
density  and  other  properties  super iaposed  on  a  one-diaensional 
flow.  This  1-D  flow  was  assuaed  to  be  achieved  by  application 
of  a  properly  selected  axial  electric  field  distribution,  which, 
when  relaxed  to  zero,  introduces  a  non-zero  Hall  current  and 
two-diaensional  effects.  The  usefulness  of  this  perturbation 
approach  is  United  to  conditions  not  too  close  to  "onset", 
but  as  the  results  contained  in  our  2nd  yearly  report  show,  the 
strong  nonuniforaities  are  likely  to  appear  only  at  currents 
upwards  of  90%  of  onset.  It  oust  be  noted  that  the  basic  flow 
upon  which  the  expansion  is  based  satisfies  the  condition  E-uB<vE, 
which  is  strictly  true  at  Ra->  *.  For  finite  Ra,  the  inequality 


applies  well  only  far  froa  the  inlet  and  exit  current  concentrations, 
as  our  later  work  has  clarified.  However,  the  convergent-divergent 
design  adopted  also  iaplies  very  weakened  inlet  and  outlet  concen¬ 
trations,  which  make  the  results  reported  in  Appendix  1  applicable 
even  at  only  moderate  Rn  values. 

The  concepts  discussed  above  are  quite  novel,  and  were  met  with 
soae  resistance.  Particular  concerns  were  raised  by  Ors.  King  and 
Lawless  about  the  validity  of  the  Alfven  (or,  more  generally, 
magnetosonic)  choking  conditions  since  it  was  their  opinion  that 
classical  sonic  choking  was  mainly  responsible  for  the  greater 
part  of  the  performance  parameters  (voltage,  in  particular)  of 
MPD  thrusters.  Dr.  Turchi  also  questioned  the  applicability  of  Rm-.>  * 
concepts  at  the  realistic  values  found  in  practice. 

We  have  investigated  the  problem  in  more  detail  since  then,  and 
by  the  close  of  the  reporting  period  had  achieved  substantial 
clarification  of  these  issues.  A  paper  is  in  preparation  for  the  1987 
Electric  Propulsion  Conference  that  will  go  in  considerable  detail 
over  the  mathesatical  structure  of  the  self-field  accelerated  plasma 
problem;  an  outline  is  included  here  as  Appendix  2.  The  essential 
result  is  that  the  channel  voltage  is,  in  fact,  largely  determined  by 
the  magnetoacoustic  .downstream)  condition  and  that  the  role  of 
classical  sonic  choking,  which  is  indeed  present  in  the  problem, 
is  limited  to  determining  the  inlet  flow  velocity  and  density. 


This  then  throws  into  question  any  explanation  of  the  onset 
phenomena  as  resulting  from  the  inability  of  the  channel  voltage 
(if  prescribed  by  upstream  conditions)  to  exceed  everywhere  the 
back  e.m.f. ,  since,  in  fact,  it  is  this  maximum  back  e.m.f.  that 
determines  the  voltage  for  given  channel  length,  current  and 
geometry.  As  to  the  applicability  of  high  Rm  concepts,  the 
answer  is  less  complete,  but  both  analysis  and  numerical  work  suggest 
that  10%  accuracy  in,  for  instance,  voltage,  can  be  obtained  with 
the  Rm-><»  limit  if  properly  used  when  Rm>10. 

The  second  of  the  tasks  for  the  year  was  the  development  of  a 
numerical  code  for  realistic  geometries  of  MPD  devices.  Progress  here 
has  been  steady  but  slower  than  anticipated,  due  to  the  difficulty  of 
the  problem.  Some  of  the  difficulty  is  clearly  highlighted  by  the 
one-~dimensional  results  just  discussed:  for  instance,  as  Appendix  2 
shows,  the  sonic  passage  condition  imposes  very  large  inlet  density 
gradients,  even  though  this  has  only  a  small  downstream  effect. 

These  gradients  will  be  very  difficult  to  resolve  in  a  2-D  grid,  but 
if  their  nature  is  understood,  we  can  hope  to  incorporate  their 
effects  in  a  manner  analogous  to  that  used  for  shocks  or  sheath  drops. 
Also  difficult  is  the  development  of  numerical  schemes  robust  enough 
to  handle  the  very  low  density  (and  therefore  high  Hall  parameter' 
regions  which  appear  near  the  anode  of  the  thruster.  These  regions 


require  careful  numerical  handling  to  avoid  appearance  of  negative 
densities  and  other  non-physical  effects  which  lead  to  computational 
instability.  More  generally,  there  is  a  general  lack  of  analytical 
guidance  on  the  stability  properties  of  the  usual  CFD  methods  when 
extended  to  our  more  complicated  situation. 

These  problems  notwithstanding,  we  have  developed  a  first- 
generation  code  documented  in  Mr.  Chanty’s  Master’s  Thesis  (Ref.  1), 
which  has  produced  meaningful  results  for  a  simple  axisymmetric 
channel  with  6  g/sec  of  Argon  at  currents  up  to  10  KA.  Convergence 
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difficulties  did  not  allow  higher  currents  to  be  included  in  this 
fully-coupled  model.  On  the  other  hand,  separate  computations  of 
flow  or  current  distributions  (with  one  or  the  other  kept  frozen.’ 
can  be  routinely  made  with  the  code.  More  recent  work  has  refined 
these  results  by  including  more  spatial  resolution,  smoothing  the 
grid  and  improving  the  iteration  method.  We  have  acquired,  with 
AFOSR  funding,  a  20*  share  in  a  multi-processor  vector  computer,  and 
plan  to  move  our  VAX  codes  to  it,  with  order-of-magnitude  increases 
in  speed  and  storage.  Appendix  3  shows  some  of  the  results  of  the 
work  in  Ref.  1. 

Turning  to  the  area  of  experimental  verification,  an  experiment 
haa  been  designed  and  built  to  test  several  predictions  of  plasma- 
dynamic  theories  which  have  been  presented  elsewhere  (see  Appendix  Ij. 
These  theoretical  results  indicate  that  area  contouring  can  be  used 
to  modify  significantly  the  performance  of  MPD  thrusters,  and,  in 
particular,  that  convergent-divergent  area  distributions  can  be 
devised  to  produce  a  nearly  uniform  current  distribution,  which  should 


result  in  reduced  frozen  losses.  Other  aspects  of  the  theory  to  be 
experimentally  tested  include  the  existence  of  magnetosonic  choking 
at  the  channel  throat  and  the  approximate  similarity  of  the  magnetized 
plasma  flow  to  that  of  an  ordinary  gas  with  a  specific  heat  ratio 
of  2. 

Two  thrusters  have  been  designed  for  these  purposes,  one  with  a 
constant  area  and  the  other  with  a  variable  area  (equal  at  the  throat 
to  that  of  the  first  thruster).  The  approximate  two-dimensional 
theory  of  Appendix  1  was  used  to  design  this  second  channel  for  a 
specified,  near-constant  current  density.  The  channels  are  show  in 
Fig.  1.  Design  conditions  are  4.25  g/sec  (Argon)  and  26  KA,  which  are 
well  below  the  onset  conditions  predicted  by  anode  starvation  theory. 
Gas  injection  is  through  four  rows  of  12  choked  holes  each,  with  six 
fast  valves  feeding  the  plenum. 

The  experimental  facility  for  the  tests  is  located  at  the 
Washington  Laboratory  of  RDA  Associates,  under  an  AFOSR  funded 
cooperative  arrangement  with  MIT.  An  existing  facility  has  been 
upgraded  over  the  last  year,  to  provide  current  pulses  of  1  msec  at  50 
KA  with  background  pressure  below  10  Torr.  Diagnostics  include 
terminal  characteristics  and  electric  and  magnetic  probe  data, 
acquired  using  high  speed  Le  Croy  digitzers,  plus  exit-plane  optical 
observations  of  plasma  density  and  velocity  and  species  temperatures 
using  a  1.26  m  Spex  spectrometer  with  0.02  A  resolution  and  both  PM 
tubes  and  a  Silicon  Intensified  Target  digitizing  camera. 


Testing  is  currently  underway  on  these  channels,  and  results 
will  be  reported  at  the  1387  Electric  Propulsion  Meeting. 

Concurrent  with  the  preceding  investigations,  we  have  under¬ 
taken  the  construction  of  a  data  base  for  MPD  technology  and  science, 
with  the  objective  of  a)  providing  easy  access  to  basic  data  and 
results  in  the  field;  .b)  creating  a  manageable  structure  onto  which 
future  results  can  be  inserted,  and  (c)  facilitating  and,  to  the 
extent  possible,  performing  cross-correlations  and  comparisons  of  data 
and  results  from  various  sources  and  viewpoints  .  We  have  so  far 
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collected  and  catalogued  a  moderately  large  number  of  publications 
and  rep'orts  (Appendix  4),  which  have  been  entered  into  a  computer¬ 
ized  file  in  D-3ase  3  ( an • Ashton-Tate  trademark),  using  an  IBM  PC-AT 
computer.  This  data  base  is  being  expanded  in  three  directions, 
namely,  by  entering  additional  reports,  by  increasing  the  number  of 
access  modes,  and  by  increasing  the  amount  of  information  per 
publication.  Currently,  the  base  can  be  entered  by  author,  date, 
laboratory  or  title.  Our  next  planned  additions  include  a  list  of 
technical  keywords  to  be  scanned,  selected  to  cover  the  MPD  field, 
and  the  inclusion  of  a  brief  abstract  with  the  paper  listing. 


At  the  same  time,  we  are  in  the  process  of  scanning  the  papers 
now  in  the  data  base  and  abstracting  our  selection  of  important 
technical  results  from  each  of  them.  This  is  clearly  the  most 
complete  form  of  the  eventual  data  base,  and  we  face  the  challenge 
of  incorporating  such  extended  information  as  tables  and  graphs  into 
the  computerized  data  base,  or,  more  likely,  of  cross-referencing 
between  it  and  a  paper-copy  supplement  containing  such  data.  In 
either  case,  the  data  must  be  classified  by  topic  rather  than  by 
source,  which  poses  an  additional  challenge.  Also,  as  we  approach 
commpletion  of  the  published  material,  we  must  dig  into  unpublished 
data  or  reports  in  the  possession  of  various  laboratories. 

3.  Personnel  Associated  with  this  Research 

The  work  was  directed  by  Associate  Professor  Manuel  Martinez- 
Sanchez,  with  the  cooperation  of  research  assistants  D.  Heimerdinger 
(Doctoral  Candidate),  J.  Marc  Chanty  (M.S.,  Sept.  1986;  currently 
a  Doctoral  Candidate),  and  Tze-Wing  Poon  (M.S.,  Sept.  1985). 
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Abstract 

A  theory  is  prsssntad  of  the  flow  of 
;■  fully  ionised  pisses  in  »n  MPD  accele¬ 
rator.  Examination  of  the  High  Magnetic 
(Reynolds  nusber  limit  reveals  the  origin 
of  the  aaode-*lip  and  cathode-end  current 
concentrations  end  shows  that  a  redistn 
button  of  current  can  be  effected  by 
contouring  the  channel.  More  generally, 
there  exists  a  distribution  of  electrode 
potentials  for  which  current  is  purely 
transverse.  Achieving  this  condition 
would  require  electrode  segment at  ion ,  but 
would  eliainate  anode  depletion  as  a  cause 
of  onset.  A  linearised  analysis  is 
presented  for  cases  with  other  electrode 
potential  distributions,  using  the  trans 
verse  case  as  a  zero’  th  order  approxima- 
tion.  The  results  show  that  strong 

transverse  aass  redistribution  can  be 
brought  about  by  even  minor  axial 

currents.  The  theory  yields  a  simplified 
onset  criterion  similar  to  that  of  Baksht, 
and  can  also  be  used  to  guide  the  design 
of  MPD  channels  with  desired  current  and 
aaaa  distributions. 
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constant  (Eq.  54) 

n  non-dimans ional ized  density 

AR  eigenvalues 

u  sheath  constant  (Eq.  35) 

ua  permeability  of  free  apace 

t  uon-d imena iona 1 i zed  mesa  flux 

*  total  preasure 

o  gas  density 

a  conduct  ivity 

r  inverse  interaction  number 

t  stream  function 

<>  transverse  average 

subscripts: 

,,  a  at  the  cathode  or  anode 

integral  index 
rtf  reference  quantity 

x,y  La  the  x  or  y  direction 

superscripts: 

0  zero’ t h  order 

’  non-dimensional  quantity:  also  used 

as  perturbation  (Eqs.  31  and  42): 
or  as  a  name  "tag"  (Kq.52) 

*  choking  condition 


Ini  roducy,on 

The  goals  of  tin*  designer  of  an  MPD 
thruster  are  to  obtain  a  desirable  overall 
performance,  plus  desirable  internal  dis 
tributions  of  electrode  current  densities, 
ohmic  dissipation,  and  possibly  gas 
density  or  asss  flux.  Purely  electrody- 
namic  considerations  can  provide  valuable 
information  for  overall  performance  (e.g., 
the  Maecker  thrust  equation),  but  the 
strong  electromechanical  coupling  must  be 
considered  in  order  to  achieve  the  goals 
related  to  internal  flow  and  currant 
distributions.  This  coupling  can  affect 
performance  in  several  important  ways  such 
as  in  the  establishment  of  the  onset 
limitation  or  m  the  determination  of 
lifetime  or  materials  choice  (from  local 
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electrode  loadings).  In  this  paper,  we 
apply  a  simple  single  fluid  aodel  of  (.he 
MPO  plasaadynaaics  to  the  problea  of 
designing  a  two-dimensional  plasaa  accel 
orator  for  prescribed  distributions  of 
flow  and  current.  The  analysis  provides 
insight  into  the  behavior  of  the  device 
and  gives  guidelines  on  Ihe  choice  and 
shapes  of  the  paraaeters  for  specific 
design  goals.  While  the  trends  and  I  mils 
obtained  are  believed  to  be  of  general 
usefulness,  their  accuracy  is  liaited  by 
the  relative  crudity  of  the  plasaa  aodel. 


Model  Formulation 

The  plasaa  is  aodelled  as  a  fully 
ionised  gas  characterized  by  a  constant 
electrical  conductivity  o,  and  a  Hall 
paraaeter  4  which  is  proportional  to  the 
aagnetic  induction  B  and  inversely  propor¬ 
tional  to  the  gas  density  p: 

3  am 

a  •  Constant  ;  4  *  c  -  (c  >  )  (1) 

P  e 

The  gas  is  assuaed  to  follow  the  ideal  gas 
law  P  -  oB(T.  with  constant  aolecular  aass 
and  teaperature  where  R«  is  calculated 
considering  an  electron-ion  fluid.  This 
last  approximation  is  partially  justifi¬ 
able  on  tha  basis  of  the  large  electron 
theraal  conductivity  and  the  large 
available  reservoir  of  ionization  and 
excitation  energy,  but  it  ia  ultiaately 
just  a  convenient  device  to  facilitate  the 
analysis. 

The  generalized  Oha's  law  is  used  in 
the  fora 

«(*  <-ux»)=j-jx4  (2) 

where  4  *  4(1/B).  The  electron  pressure 
gradient  tera  is  neglected  because  Baksht 
et  al.  (Reference  1)  included  this  effect 
in  their  aiaplified  analysis  and  found 
that  the  results  were  siniaally  aodified, 
with  no  effect  on  the  overall  behavior. 
The  ion  slip  is  also  neglected,  due  to  the 
fact  that  the  electron  Hall  paraaeter 
rarely  exceeds  values  of  the  order  of 
unity  io  MPO  channels. 

Continuity,  aoaentua,  and  the 
pertinent  Maxwell's  equations  are  listed 
below: 


4^pu )  #  4^pv )  .  Q 

(3) 

3  * 

*(«•*)«  »  -  V  (  p  *■  * - ) 
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(4) 

j  -  y  x  !-  ;  Bill) 

(5) 

Me 

V  x  1  -  0 

(6) 

and  j  can  be  eliainated 

and  an 

induction  equation  foraed 

4*8  ^  4*8  f  4B  44  _  4B  44  ] 

4x*  4 yd  l  Jy  ~  7x  4y  J 

s  (  |  ]  -  p»oo  •  (  uB  )  (7) 


where  the  Laplacian  teras  represent  I  lie 
aagnetic  diffusion,  the  third  tera  arises 
froa  the  variations  of  the  Hall  paraaeter 
due  to  its  inverse  dependence  on  gas  den¬ 
sity,  and  the  right  hand  side  represents 
the  effect  due  to  aagnetic  convection. 

It  is  useful  to  non-d laens ions t i ze 
these  equations  using  reference  quantities 
representative  of  those  encountered  in  MPO 
channel  flow.  taposition  of  continuity 
and  aoaentua-  like  constraints 
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lOrefUrefh*  -  (j  and 
.providas  for  these  aeunmgful  reference 
paraaeters.  This  results  in  the  non- 
diaenaional  variables 


- .  v 

,  y 

.  vL 

*  •  r 

>  Fa 

V  '  — '■ 

Urn 

u.  *  j i _ 

««  ,  £_ 

B‘  =  b 

Uref 

Or  *4 

.8) 


and  tha  non-diaensional  equal  ions 
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The  following  diaansionless  paraaeters 
appear: 
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The  High  Magnetic  Hevnolda  Number  Liai l 


Even  though  R«  is  not  likely  to  be 
really  large  in  MPD  channels,  it  is  large 
enough  that  some  information  on  the  flow 
behavior  can  be  gained  by  considering 
first  the  limit  Re  >>  1. 

Froa  this  point  on  we  will  drop  the 
(')  to  denote  diaansionless  quantities. 
As  R*  w  •  and  v  -  0,  B/p  becoaes  a 
convective  constant.  The  effective  total 
prtsiur*  is  now 


rp 


if  i  )'- 


(  17  ) 


and  is  a  function  of  density  alone.  Being 
barotropic  aeans  that  a  Bernoul l i-type 
conservation  equation  can  be  obtained: 


3MI)l)WMMt)«i)l)^^^ 
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(•mini  current  is  ssoothly  distributed 
along  the  nozzle  with  current  lines 
corresponding  to  constant-0  lines. 


«••  (  3*  ►  r  m  .  .  [  |  j*.  ]  - 

showing  that  the  quantity  in  brackets  is 
another  convective  constant. 

tn  this  liait,  the  general  behavior 
>*f  the  plans*  flow  can  therefore  be 
described  as  being  a  direct  analog  of  an 
ordinary  coapressible  flow  with  the 
particular  pressure-density  variation  ». 
which,  if  r  ■.  I  is  nearly  that  f or  a 
spot  i lie  heal  ratio  v  ■  2,  If  the 
convert i ve  constants  are  constant  across 
»t  renal  tnes  then  all  the  taailiai 
properties  of  two -diaens tone l  ideal  gas 
flows  then  have  their  counter-part  as  high 
Magnetic  Reynolds  nuaber  I  lows:  choking, 
I’ranill  l -Ms.ver  expansions,  shocks.  sir. 
The  role  of  the  speed  of  sound  is  here 
ployed  by  the  disturbance  speed 


which  is  the  aagneto -acoust i c  velocity  and 
reduces  to  the  speed  of  sound  at  zero  Bag- 
net  ic  field  and  to  the  Alfven  speed  for 
low  pressure  and  high  field. 

Two  laportent  identifications  can  be 
made  in  these  flows: 

(a)  Iso  density  lines  are  also 
iso  field  lines,  hence  they 
are  electric  current  lines, 
lb)  Straaalinea  are  also  equi 
potentials. 

Consider  a  crude  aodel  of  an  MPI)  chan 
nel  consisting  of  a  rectangular  channel 
discharging  into  a  vacuus  where  R  =0  (Fig. 
1).  Since  the  cons  tan t -area  section  dis 
charges  into  a  vacuus,  it  sust  be  a  region 
of  critical  flow  lu  -  si.  Pros  the 
Bernoulli  integral  iHq.  18),  the  sagneto- 
icouitic  speed  (Eq.  19),  and  the  asausp- 
tions  r<<l  and  u(0)«0,  u= s- l ( 2/3 ) /  p% ) ‘ '* 
which  is  f  2 / 3  )  1  ^  *  tines  the  stagnation 
Alfven  velocity.  Also  8  •*  2/3  (B(0)  i  li 

and  o  =  f2/3)0i 0) .  The  flow  suet  sake  an 
abrupt  transition  f ros  p,  to  (2/3)s«  as  it 
accelerates  rapidly  to  the  critical  speed 
and  the  corresponding  fall  of  R  isplies  a 
concentrated  current  sheet. 

The  expansion  into  vacuus  proceeds  by 
way  of  a  pattern  of  lip-centered  expansion 
tans  which  interact  eutuallv  downstream  to 
fora  a  supercritical  jet.  The  expansion 
fans  also  represent  current  concentra¬ 
tions.  Fully  2/3  of  the  total  current  is 
carried  by  these  pseudo  arcs,  and  the 
current  pattern  in  the  jet  can  be  predic¬ 
ted  by  the  well  established  aethod  of 
characteristics.  For  r  ••  <  1,  a  siaple 
calculation  gives  for  the  liaiting  sirens  • 
line  land  also  the  liaiting  Mach-Alfven 
line)  an  inclination  of  ( V'J  -  11/2 
radians,  or  63.8*  (Fig.  2). 

If  a  convergent - d i vergent  nozzle  were 
provided,  the  "chamber"  flow  would  be 
subcritical,  since  the  critical,  or 
eagneto-acoust ic  speed  is  reached  exactly 
at  the  throat  and  the  exit  flow  will  now 
be  super-cr i t ica 1 .  The  initial  juap  in  B 
is  now  decreased  and  the  exit  iip  current 
concentration  it  also  weakened.  The  re- 


The  Slender  Channel  Approximation 

Me  wish  to  exasine  the  flow  in  cases 
with  sore  realistic  paraseters,  with 
esphasis  on  providing  a  asooth  axpanaion 
that  say  yield  a  sore  benign  current  die 
tnbutinn.  Me  will  exploit  the  slender- 
nose  of  the  likely  geometries  by  assusing 
that  h*/L  <  '  l . 

It  in  conveniant  to  usa  as  the  trans¬ 
verse  independent  variable  the  stress 
funct  ion  t,  defined  bv 

9t  d*  , ,n 

ou  '  3y  0v  -  -  ^  .20, 

where  t  is  nond leans  tone l czed  by  G.  The 
governing  aquations  reduce  to 
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where  3  R*  ^  and  3  -  8r#i  ^  s l 24 ) 

To  ispoee  the  electrode  well  boundary 
conditions,  Oha’s  law  is  solved  for  the 
axial  field,  E*.  (normalized  by 
resulting  in 
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(25) 


The  implication  of  Eq.  (25)  is  that 
the  I ransverae  Lorentz  forces  due  to  any 
axial  current  are  hydrostatically  balanced 
by  a  lateral  pressure  gradient. 


Design  ./or.  Zero.  Axial  Currant 

Thera  is  interest  in  eliminating  the 
axial  current  due  to  the  Hall  effect  which 
has  been  related  to  anode  starvation  and 
onset  Ref  l).  Examination  of  Eqa. 
(21)-(23>  shows  that  solutions  do  exist  m 
which  all  the  variables  depend  on  x  alone. 
Again,  two  convective  constants  are  found 
and  explicit  solutions  can  be  found  given 
initial  values  of  p,  B,  and  an  axial 
pressure  specification  of  x(x).  These 
expressions  are 

»  2  iB/pi1*,*)  -  r 

0  — - 

8/0  )  * 

l  '■  S7ST  ;  8  =  (  l  ]  0 

U  --  |  u  l  0  )  1  *  2  r  In  ~  2  (  B/0  )  1  ( 0  (  0  )  -0  ) 
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At  the  "choking  point",  the  locel 
density  p *  is  given  by  the  solution  to  the 
equation  u  *  s,  or 
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The  height  distribution  is  shown  in 
Fig.  4  as  s  function  of  8/S(0)  for  r  •  0. 
For  the  particular  case  when  a  constant 
current  density  is  desired,  the  curve 
shows  the  true'  shape  of  the  channel. 

The  price  that  eust  be  paid  for  the 
achieveaent  of  a  purely  transverse  current 
is  the  required  presence  of  an  axial 


electric  field. 
6q.  (25): 


Its  value  follows  from 


This  isplies  the  iaposition  of  an  axial 
field  through  a  segaented  channel  with 
independent  control  of  the  voltage  of  each 

segaent . 


Current 


L-_  Foraulat  ton 


When  the  axial  field  departs  froa  the 


value  given  by  Bq. 


axial  currents 


and,  consequently,  lateral  density  gradi- 


ants,  will  appear  in  the  plaaaa.  A 
particular  exaaple  is  that  of  a  channel 
with  continuous  aetallic  walla;  if  the 
axial,  variations  of  the  electrode  voltage 
drops  can  be  neglected  in  coaparison  to 
the  Hall  field  of  Bq.  (30),  then  E«  will 
be  zero  along  the  walls,  and  j*  will  be 
nonzero.  When  the  voltage  drop  variations 
are  accounted  for,  Ea  will  not  be  precise¬ 
ly  zero,  but  will  not  generally  agree  with 
Eq.  (30)  either. 

We  will  assuae  relatively  saall 
departures  froa  the  purely  transverse 
condition  where  BK  is  given  by  Eq.  (30). 
The  quantities  corresponding  to  this 
zero'th  order  solution  will  be  denoted  by 
a  zero  superscript,  and  the  equations  will 
be  linearized  in  the  nondiaens 1  one  1 
perturbation  quantities,  defined  by 


a  -  p®(l  *  n)  i  r  -  0 u  =  r®( 1 

0  *  !•  *  B* 


but  with  no  perturbation  in  *.  In  addi¬ 
tion,  the  channel  profile  is  to  be 
perturbed  to  accoaodate  the  0 u  pertur¬ 
bations  brought  on  by  the  Sx  variations. 

Performing  the  linearization  of  Sqs. 
(21)and  (25)  gives  new  perturbation 
induction  and  electric  field  equations 
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The  coefficients  m  Bq.  (32 1  and  i33< 
are  functions  of  x  alone,  and,  m  fact, 
B°/p®  is  a  pure  constant  equivalent  to 
(B/p)®.  The  factor  (dw/dx)/B®  would  be 
nearly  constant  in  a  cons  tan t -cur ran t 
zero’th  order  channel,  and  r®  is  flat  in 
the  throat  region  becoaing  saall  as  the 
channel  diverges.  Thus  we  are  justified 
in  adopting,  for  an  approxiaate  solution, 
the  axially  averaged  values  of  the  coef¬ 
ficients;  Eq.  (32)  will  be  written  as 


where  a  and  b  are  constants  given  by 


:  b  .  3  ,U  -  g  (35) 

H(  B/j# )  •  *  H  1  1  Jx  J 


the  synbol  <  >  indicates  axial  averaging. 

For  analysis,  we  will  assuae  that  the 
perturbation  K«*  of  B«  on  the  electrode 
walls  is  a  prescribed  fraction  t -p >  of  the 
zero'th  order  value,  so  that. 
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where 


<  r®  > 


3  ( B/p ) • 


(B/p)®  .  dB® 

e  ,  - -  <-  . 


The  precise  value  of  p  eust  be 
obtained  fros  a  nodal  of  the  sheaths  on 
the  electrode  walls.  If  we  assuae  metal- 
lic  walls  with  constant  voltage  drops  so 
that  the  total  E*  is  zero  at  the  walls,  »e 
eust  use  u  -  *1.  This  is  the  only  rase  we 
have  so  far  analyzed. 

In  adopting  the  slender  channel 
assumption,  all  diffusion  terns  were 
neglected  Irnnsforning  the  problem  from 
elliptic  to  parabolic.  Therefore  an  extra 
condition  on  n  <o(x-0,p))  nust  be  pre¬ 
scribed  but  no  condition  can  be  inposed  on 
downstreaa  boundary  conditions. 

In  a  aatheaalical  sense,  we  are  sol¬ 
ving  only  the  outer  liait  of  a  singular 
perturbation  problem  with  the  saall 
parameter  (h*/L)*.  The  inner  outer 
astrhmg  requires  the  eliaination  of  the 
exponentially  growing  parts  of  the  outer 
solution.  We  have  investigated  the 
question  of  whether  a  similar  inner-outer 
Batching  is  required  near  the  x  -  0  end, 
i.e.,  whether  our  "outer”  solution  can  be 
aade  to  fit  the  initial  conditions  or,  on 
the  contrary,  a  near-d i scon t 1 nu i t y  may 
appear  at  x  =  0.  The  results  of  this 
analysis  indicate  that  no  "boundary  layer” 
is  forced  at  x  x  0  on  the  zero'th  order 
solution  as  long  as  no  axial  current  is 
allowed.  Otherwise  the  diffusion  tera. 


d*/r/dx*.  cannot  be  neglected.  Notice  that 
the  abaanca  of  a  strong  tranaition  at 
x  -  0  is  not  a  universal  raault;  in  fact, 
it  we*  previously  seen  that  ona  such 
singular  layer  does  appear  when  the  araa 
distribution,  i os  toad  of  tha  total 
praasura  distribution,  ia  prescribed. 

To  coapiele  tha  linearized 
foraulation,  tha  perturbation  on  tha  aass 
flux,  £<|.  (22)  ia  calculated: 


d  In  u* 


which  integrates  to 


in  d  In  ua 


t(x,T)  -  flO.P)  <•  n  -  /J(O.T)  - 

“  U#T  <o  "  dx  t40) 

with  I  -  <  -  i,  >  B  (  5  j%  141) 

Tha  channel  height  follows  froa  tha  inte¬ 
gration  of  tha  straaa  function;  the 
zero' th  order  value  ia  h*  *  l/f*,  and  ita 
perturbation  is  given  by 

srfxT  *  "  *  d*  ~  *  rnm  *  a  - 


-  nTVTTT  -  jfr  5  dx 


where  the  overbar  indicates  trana.verae 

averaging. 


2-— -SalJLt.ian_  of  the  .  Linearized  gquetions 

Th*  solution  to  th*  above  aquation* 
for  th*  cas*  of  equipotential  walla  (x*\ 
in  Bqi  (36))  is  presented.  Using  the 
nethod  of  separation  of  variables,  and 
taking  special  car*  to  acconodate  the 
presence  of  axial  derivatives  in/ix  in  the 
boundary  conditions  (Kq.  (36)),  w*  arrive 
at 

n  -  l  *  C. i exp( -A. ix) exp( r. i t)  * 

b  • 

*  expt-  T)  T  C„exp( -A,x ) ( cos  n»f  - 
-  Kaain  nxT)  (43) 


an  additional  C.j  tern)  go  to  zero.  Th* 
sinplest  type  of  initial  density  profile, 
with  zero  aean,  capable  of  satisfying  this 
condition,  is  linear 

h(O.T)  -  n<;(  2*  l  )  and  (47) 
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Th*  coefficients  C.|,  CB>«  aust  be 
obtained  by  integration  perforned  using  an 
initial  density  distribution  n(x-0,t),  and 
by  satisfying  a  nodified  Fourier  orthogo¬ 
nality  condition.  Th*  initial  density 
distribution  aust  also  be  chosen  so  that 
all  axially  divergent  terns,  (in  this  css* 


3 j _ A  Nunerical  Exasol* 

Th*  following  exanple  was  chosen 
aainly  for  purposes  of  illustration  and 
for  verification  of  th*  theory.  Its 
paraneters  bear  son*  resenblanc*  to  those 
encountered  in  Argon  MPD  thrusters  opera¬ 
ting  near  onset  with  several  grass  per 
second  of  flow.  We  lake 

R,T  =  2.496x10*  J/kg  B • ( 0 )  *  0.05  Tesla 
c  .  0.00207  kg/a*/Tesla  a  -  5000  aho/a 


u*(Q)  *  3000  a/s  0 

«*(0)  1 . 333*10**  kg/a* 


0.02  kg/aec/a 
I.  0 .  I  a 


-  »««>(  1  -  J7T  Y 


with  *  is  aatars. 

Tha  solution  for  saro  axial  currant 
is  presented  ia  figures  3  and  6  in 
diaaaaioaal  taraa. 

for  tha  purposas  of  tha  linaartzad 
solution,  wa  adopt  approxiaata  averages  as 
foil owo : 


<r>  *  0.4  kg/sec/a* 


<2^  =  TO  V/i 


<  -  ~  >  -  2.090  x  10*  a/sac* 


With  thasa  valuaa,  tha  linaarized  aqua 
tiona  giva  n  expansion  coefficients  (Table 
1)  and  tha  sala<-t«d  results  shown  tn 
figuraa  T,  8.  and  9.  figure  7  shows  that 
tha  total  dansity  p-p*(l >n)  is  strongly 
concantratad  near  tha  cathode.  to  the 
point  that  slightly  negative  densities  are 
calculated  near  tha  anode,  particularly  in 
tha  region  near  and  just  beyond  tha 
throat.  This  is  a  clear  indication  that 
tha  channel  has  reached  onset  conditions. 
Tha  aass  flux  au  has  bean  assuaed  uni  fora 
at  x  *  0,  and  raaains  not  fnr  froa  uni  for 
aity,  which  forces  tha  velocity  to  be 
substantially  higher  near  tha  anode. 

tn  figure  8  wa  sea  the  relatively 
ainor  transverse  nonun i forai ty  of  ft  which 
ia  raaponaibla  for  the  large  nonuniforai- 
tiea  of  a  and  u  shown  in  figure  7.  The 
Mall  paraaater  (which  is  approxiaatel v  the 
ratio  of  ja  to  j*.  since  EK  is  now  near 
taro  everywhere)  is  also  shown  in  this 
figure,  and  is  large,  near  tha  anode. 

figure  9  coa pares  the  unperturbed  and 
perturbed  channel  contours.  Tha  differen¬ 
ces  ara  aodarata,  indicating  that  the 
average  aaas  flux  is  only  weakly  effected 
by  tha  aasa  redistribution. 

As  shown  in  figure  8,  8  bacoaaa  zero 
on  tha  cathode  at  about  x  =  0.079  a,  while 
B  is  still  positive  at  tha  anode 
(B*  ^  0.017  Tesla).  Aasuaing  this  is  the 
actual  channel  taraination  point,  the 
raaaining  ■  at  the  anode  aust  drop  to  zero 
through  a  concentrated  anode  lip  "quasi- 
arc”,  which  carries  about  34k  of  the  total 
current.  for  coaparison,  truncating  the 
channel  at  x  *  0.079  a  on  the  design  with 
zero  axial  current  would  hai(e  left  anode 
and  cathode  lip  currents  of  about  I7h  of 
tha  total.  Clearly,  the  effect  of  the 
Mall  paraaater  has  bean  to  strengthen  the 
anode  lip  concentration  while  unloading 
tha  cathode  tip. 


Wa  have  indicated  that  soae  looseness 
raaains  in  the  specification  of  the  ini¬ 
tial  density  profile  for  coapa¬ 
rison  with  the  previous  exaaple,  where  a 
zero-aaan,  linear  function  was  used,  we 
use  in  this  section  an  even  ampler 
approach,  which  will  also  yield  a  ampli¬ 
fied  criterion  for  onset.  To  this  end,  we 


return  to  Bu.  (43)  end  select  n(0,P)  such 
as  to  be  orthogonal  ( in  the  aodified  sense 
required  by  tha  specinl  boundary  condi¬ 
tions  of  our  problea)  to  alt  the  nodes  for 
n  k  l.  This  leaves 


n  -  l  ►  C„ | exp<  A. i x )exp( r. i t)  (56) 


and  for  this  to  have  zero  aean  value 


expi r . i )  - 


Using  tha  above  Until,  C-i  -2.983,  coa 
pared  to  -  2.661  for  the  linear  profile. 

Tha  density  at  the  nnode  root  ix  () , 
*  -  0)  ia  given  by 


2  ♦  C.i 


Equating  this  ml io  tu  zero,  is  a  staple 
criterion  for  onset  since  this  is  the 
first  spot  where  the  density  becoaas  zero. 
Therefore  C.|*-2,  or  r.|--l.593.  The 
quantity  r. ,  is  still  too  coaplex  to  give 
any  real  inaight.  But  to  a  crude  first 
upproxiaation,  tha  (b-r)  liras  in  Cq.  t 44 i 
can  be  neglected,  leaving  r.(  »  (e/a')*'*. 
Accepting  this  gives  the  onset  criterion 
e/s  -  2.94. 

froa  Rqs.  (39)  and  (37),  (since  the 
nond taans i ona 1  variables,  like  B*,  a*,  f*, 
etc.  nr*  all  of  order  unity) 


Ms  h* 

77  r 


Using  the  xpprox mat  ion  *  waI/W 

and  nlso  C  -  »  w,  we  obtain 


—  M*'*  *  1.2 


w  !.*'*•  WT'*'* 


'  fO*  ' 1 Us  1 


Equation  (60)  displays  aost  of  the 
observed  trends  is,  M,  length,  etc...). 
Aside  froa  the  nuaerical  accuracy,  this 
expression  contains  the  suae  basic  physics 
as  that  derived  in  Reference  (li. 


Cone  1  us i on 


It  has  been  shown  I  hut  appropriate 
taylonng  of  the  channel  shape  in  an  MP0 
accelerator  can  be  used  to  obtain  a 
desired  axial  distribution  of  transverse 
current,  and  that  taylonng  of  the  elec 
trode  potential  distributions  can  be  used 
to  control  and  suppress  axial  currents  xnd 
lateral  density  gradients.  While  the 
practical  iap I eacn t o t i on  of  these  notions 
say  present  difficulties,  the  theory 
presented  can  be  used  as  a  guideline  for 
designs  with  perhaps  better  perforaance 
and  longer  lifetiaes. 
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a 


A 


a 

.  A« 

<• 

c. 

•  2 

-  24.332 

0 

l 

0.4450 

2.664 

l 

4.441 

0.8683 

0.9070 

2 

14.992 

1.0834 

0.0371 

3 

32.470 

1.4437 

0.0370 

4 

52.996 

1.8402 

0.00714 

S 

68.530 

2.2513 

0.00627 

6 

122.07 

2.6696 

0.00184 

7 

172. 6H 

3.0920 

0.00180 

H 

225.  17 

3.5171 

0.00064 

9 

284.7 

3.9439 

0.00069 

10 

351.3 

4.3719 

0.00028 

11 

424.9 

4.8008 

0.00032 

12 

505.5 

5.2303 

0.00014 

13 

593.  1 

5.6603 

0.00016 

14 

687.7 

6.0908 

0.00008 

19 

789.3 

6.5315 

0.00009 

Table  1:  Coefficients  for  Eq.  43 
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Fig.  1:  Flow  it  la  «  •  Froa 
Rectangular  Channel 
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fig.  3:  Saooth  Expansion  and  Current 
Distribution  in  Convergent  - 
Divergent  Nossie 


Fig.  4:  Nozzle  Contour  for  r  *  0 


Fig.  2:  Corner  Expansions  Into  Vacuus 

a.  Critical  Approach  Speed 

b.  Supercritical  Approach  Speed 


Fig.  5:  Channel  Height,  Density,  sod 
Velocity  Profiles  for  Zaro 
Axial  Current 
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fit-  6:  Magnet ic  Field.  Transverse 
Current,  and  Axial  Field 
Profiles  For  Zero  Axial 
Current. 
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Fig.  7:  Transverse  Oensity.  Velocity  and 
Mass  Flux  Profiles  at  Several 
Axial  Oistances,  For  Channel 
with  Iquipotent tal  Walls 
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Fig.  8:  Transverse  Profiles  of  Magnetic 
Field  and  Hall  I'araaeter  at 
Various  Axial  Distances,  for 
Channel  with  Equ  i  po t en 1 1 a  1  Walls 


Fig.  9:  Channel  Height  Profiles  for 
Unperturbed  i j,  .01  and 
Perturbed  0)  Cases 
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THE  STHUCTURB  OF  SELF-FIELD  ACCELERATED 
PLASMA  FLOWS 


By 

Manuel  Mart inez-Sanchez , 

MIT,  Departaent  of  Aeronautics  L  Astronautics 


This  paper  will  present  a  detailed  analysis  of  the  developaent  of 
flow  and  electroaagnetic  properties  in  a  contoured  quaa i-one-diaen- 
sional  MPD  channel,  with  a  clarification  of  the  roles  of  ordinary  sonic 
passage  and  aagtie t oacous tic  choking,  and  an  assessaent  of  the  effects 
and  potential  uses  of  area  variation. 

The  tipproach  starts  with  the  foraulation  of  a  set  of  governing 
equations  incorporating  the  condition  of  zero  axial  current,  which  is 
required  for  a  quaa i-one-diaens ional  aodel  to  be  consistent,  and  which 
iaplies  a  certain  aaount  of  axial  electric  field  (although  the  results 
are  also  representative  for  axially  shorted  channels  at  conditions 
reaoved  froa  the  starved  reg'iae).  Let  a  be  the  mass  flow  rate,  A*  the 
throat  area,  L  the  channel  length  and  Bg  the  initial  magnetic  field, 
and  define  nondiaensional  density,  pressure,  enthalpy,  velocity, 
magnetic  and  electric  fields,  current,  area,  and  distance,  by 


P  - 


h 


( ±£&J5)t 
A*  B* 


h 


19 


u  =  «_  u 


A*  BJ 


-  o_  2^ 

b  am  ai  a» 


A'  B* 


»’h  J*t*J 


A  .  _  x 

•  *  IF-  «  *  L 


Than  we  obtain 


p  u  a  =  1 


i5^*ac(5*b')  =  0 


2  “f  *  s3«  '  2*J,;  s  = 


aJ  =  E  -  ub 


—  d  In  a 

"b  ~dT~ 


Where  the  paraaeter  a  is  the  inverse  of  the  magnetic  Reynolds  nuaber 


1  _  oLA  Bii 

a  *  2a 


Nuaerically,  *  ~  0.1  -  0.3  for  realistic  MPD  thrusters.  For  small  a, 


the  above  systaa  becomes  singular,  indicating  the  potential  for  initial 
and  final  "layers'*  with  rapid  change  of  the  properties  (including 
current  concentration).  A  formal  expansion  in  powers  of  a  leads  to  an 
"outer"  set  of  equations  where  (3)  gets  replaced  by 

*®Q®  ^Si  -  5o  4|L  =  o  (8) 

i.e.,  isentropic  flow,  and  (4),  (5)  and  (6)  simplify  to 

^  (u®b® a)  =  0  (9) 

which,  together  with  (1),  indicates  b® /p®  =  const..  This  "outer"  flow 
can  then  be  reduced  easily  to  ordinary  nozzle  flow  with  an  equivalent 
"t"  value  close  to  2,  but  modified  slightly  by  the  p  term  in  (2). 

On  the  other  hand,  regions  near  the  channel  inlet  and  exit,  where 
specific  boundary  conditions  have  to  be  met,  constitute  the  "inner" 
regions,  with  thickness  4(  ~  a.  In  these  regions  we  can  ignore  area 
variation  and  obtain  the  classical  integrals  (e.g.  King,  19tl . ) 


E(  *  §o  =  const.  (10) 

Pi  u,  a<j  =  1  (11) 

u,  ♦  ao  (p;  ♦  bf )  =  u,  *■  ao  (po  +  1)  (12) 

\  ♦  |  u?  +  EaoEob;  =  hr  «■  2a<,Eo  (13) 

&  o 


together  with 
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1 


f  *  -  (i  '  5->» 


Analysis  of  this  inner  set  of  equations  shows  that  sonic  passage, 
if  it  occurs,  requires  a  special  saoothness  conditions,  which  for 

Y  •  1  • 

constant  t  is  u,b,  =  ■—  ■  ■  Ea .  In  addition,  the  requirement  that  inner 
and  outer  regions  merge  asymptotically  gives  a  specific  outer  boundary 
condition  for  the  inner  variables,  i.e.,  bj,  =  Eo  (where  "1" 

g  — 

represents  the  overlap  region).  Given  a,  A  ,  Bo  ,  hr  and  L,  a  count  of 

0 

equations  and  unknown  shows  that  in  order  to  complete  the  initial 
region  specification,  a  downstream  boundary  condition  on  the  outer 
solution  is  also  needed.  This  is  furnished  by  the  condition  (which  can 
be  inferred  from  the  outer  set,  Eqs.  (1)  -  (6>;  of  magnetoacoustic  flow 
at  the  throat  (or  throughout  the  length  if  it  is  a  constant-area  duct): 


where  c*  is  the  sonic  speed  (J r  ^  in  a  constant-^  model). 

*  P 

On  the  other  hand,  if  a  condition  b  =  0  is  imposed  at  the  channel 
exit,  as  a  crude  way  of  "bottling  up”  the  fringe  magnetic  field,  then 
no  extra  conditions  need  to  be  imposed  on  the  outer  solution,  other 
than  its  smooth  matching  to  the  new  inner  solution,  which  leads  to  b  = 
0  through  a  second  -  «  layer. 

These  equations,  when  solved  together,  give  the  first  order  (in  *) 
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solution  for  the  remaining  unspecified  inlet  condition  (pq  ,  for 
instance),  the  electric  field  Eo  (and  E($)  sore  generally),  and  other 
quantities  of  interest  at,  for  instance,  the  sonic  point,  the  end  of 
the  inlet  high  current  layer,  the  magnetosonic  throat,  or  the  channel 
exit.  We  have  obtained  solutions  for  both,  a  constant— r  model,  and  an 
Argon  ionization  equilibrium  model  of  the  gas,  and  verified  the  two 
models  against  each  other  in  low  temperature  cases.  A  sample  of  these 
results  is  shown  in  Table  1,  corresponding  to  the  specification 

■H 

=  0.50  kg/sec/m*  B0  =0.1  Tesla 

A 


Tr  o  t  .  a  = 

400  K 

A* 

=  2 

Location 

1  b 

1 

i 

u 

3 

i  ? 

.! 

T3 

.  Mach  # 

i 

1 

Mach- 
Alfven  # 

_ 

Inlet 

j  1 

0.0143 

34.882 

! 0.0443 

1 

0.2832 

1  0.3114 

0.0588 

Sonic  Point 

\  0.9793 

0. 4.157 

'”473225 

,  0.0347 

0.2832 

1 

0.1711 

Edge  of  inlet 
layer 

0.9401 

i 

1.6599 

0.2832 

1 

2.2987 

0 . 2896 

Throat 

j  0.6527 ;  0.8677 

1 

1.1525 

0.0093 

0.5664 

7.4775 

1 

Exit 

(Infinite 

Expansion) 

i 

1  0 

1.5073 

0 

i 

i  o 

1 _ 

0 

' 

GD 

_ J 

• 

_  _  '  _ i _ i _ 1 _ ! _ j._. 

TABLE  1.  NONDIMENS IONAL  QUANTITIES  FOR  CONV-DIV  MPD  CHANNEL 


Several  features  of  Table  1  are  noteworthy: 
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a)  The  flow  is  sub-magnetosonic  (but  mostly  supersonic)  in  the 
"chamber”  part  of  the  channel. 

b)  There  is  a  weak  current  concentration  at  inlet  (carrying  1  - 
0.9401,  or  6%,  of  the  total  current). 

c)  Because  of  this  last  feature,  the  actual  Argon  temperature  in  the 
chamber  is  only  3100  K,  too  low  for  equilibrium  ionization. 
Non-*equilibrium  effects  should  be  accounted  for,  but  this  result 
emphasizes  the  "off-loading"  of  the  channel  inlet  area. 

d)  There  isi drastic  density  drop  through  even  this  weak  current  layer. 
This  will  hamper  efforts  at  straight  numerical  simulation. 

e)  The  ratio  p  of  pressure  to  magnetic  pressure  remains  at  a  few 
percent  or  below  throughout  the  flow. 

For  comparison,  Table  2  shows  the  rusults  for  a  non-con toured, 
constant-area  channel  wiiV.  otherwise  the  same  data.  The  entire  channel 
(except  for  the  inlet  and  exit  "layers")  is  now  magnetosonic.  The 
inlet  region  carries  (1  -  0.63),  or  37\  of  the  current,  the  rest  being 
carried  by  the  exit  current  concentration.  The  extra  inlet  dissipation 
raises  the  Argon  channel  temperature  to  7660  K  (2.3S  ionization, 
Ttff  =  1.18)  but  it  (and  the  exit  concentration)  lowers  the  efficiency, 
as  shown  in  the  exit  nondimens ional  velocity  of  about  1.12  (1  from  EM 
effects,  0.12  from  electrothermal  expansion),  compared  to  1.51  for  the 
case  in  Table  1.  In  fact  for  the  perfect  expansion  of  Table  1,  the 
thermal  energy  resulting  from  the  inlet  dissipation  is  fully  recovered, 
and  the  thrust  effeciency  is  unity;  but  even  if  we  assume  no  recovery 
beyond  the  throat,  an  efficiency  of  0.95  is  calculated  (ignoring,  of 


course,  all  transport  losses).  By  comparison,  the  thrust  efficiency  of 


the  constant-area  channel  is  only  0.687.  Note  also  the  even  more 
extreme  density  drop  through  the  inlet  layer  in  this  case. 


Location  b  u  p  p  E  Mach  #  Mach/ 

Alfven  * 


Inlet 

HBi 

ISSl 

0.1520 

0.4576 

07187 

0.0622 

Sonic  Point 

0.9304 

0.1786 

5.598 

0.1156 

0.4576 

1 

0.3058 

Channel 

0.6274 

6.7290 

1.3717 

0.360 

0.4576 

4.188 

1 

Exit  Plane 

(B  =  0) 

0 

1.1212 

0.8919 

0.0374 

0.4576 

5.238 

5.238 

TABLE  2.  NONDIMENS IONAL  QUANTITIES  FOR  CONSTANT  AREA  CHANNEL 


More  complete  results  and  trends  will  be  included  in  the  full 
paper,  including  extension  to  higher  order  in  inverse  Magnetic  Reynolds 
number  and  comparisons  to  straight  numerical  simulations.  Work  is  also 
underway  on  the  related  issues  of  wall  losses  and  non-equilibrium 
ionization  effects.  Two  companion  papers  will  present  two-dimensional 
simulation  effects  obtained  numerically,  and  experimental  results  on  a 
pair  (contoured-straight )  of  channels. 
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Two-Dimensional  Numerical  Simulation  of  MPD  Flows 


by 

J.  Marc  Chanty,  S.C.  Myer  and  M.  Mart inez-Sanchez 

MIT 


This  paper  will  present  results  of  our  on-going  work  on  numerical 
aodeling  in  two  or  three  dimensions  of  the  fluid  and  electrodynamic 
property  distributions  that  occur  in  plasma  accelerators.  This 
numerical  approach  is  an  important  supplement  to  experimental  and 
analytical  studies,  both  of  which  have  limitations  of  their  own;  yet, 
application  of  numerical  models  to  electric  propulsion  has  lagged 
behind  other  fields,  where  such  applications  are  now  more  or  less 
routine.  In  part  this  is  due  to  the  inherent  difficulty  of  the 
governing  equations  where  fluid  and  electromagnetic  fields  must  be 
allowed  to  interact  strongly,  with  the  result  that  several  spatial  and 
temporal  scales  occur,  convective  and  diffusive  effects  of  several 
types  compete,  and  the  classical  subsonic-transonic-supersonic  regimes 
of  gasdynamics  must  be  expanded  due  to  the  occurrence  of  characteristic 
wave  speeds  other  than  that  of  sound. 

We  have  developed  a  two-dimensional  axially  symmetric  model  of  an 
MPD  thruster  using  a  finite-volume  discretization  approach  in  order  to 
accomodate  the  complex  geometry,  an  example  of  which  is  shown  in  Fig. 
1.  An  automatic  grid-generation  code  has  also  been  produced  to  help  in 
this  process.  The  conservation  equations  for  a  1-fluid  model  of  the 
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plasma,  and  the  Maxwell's  aquation*  are  satisfied  in  the  steady  state 
by  balancing  the  appropriate  fltacee  into  and  out  of  each  of  the  finite 
voluaee.  The  equations  are  then  linearized  successively  about  the  beat 
currant  solution,  and  incremental  corrections  are  calculated  throughout 
the  flow  by  solving  e  large  set  of  linear  equations,  the  coefficients 
for  which  fora  a  banded-structure  matrix  that  is  updated  after  each 
linear  solution.  This  Newton-Baphson  iteration  converges  rapidly  (less 
than  10  cycles),  but  requires  a  large  amount  of  data  storage  (about  3 
Mb  for  a  20  x  40  grid). 

Stability  of  the  numerical  process  requires  satisfaction  of 
several  criteria  similar  to  those  known  from  ordinary  CFD  work,  but  a 
complete  study  of  the  stsbility  problem  for  our  more  complicated  case 
ia  still  lacking.  We  find  empirically  that  artificial  damping  terns 
must  be  introduced  in  the  equations,  the  effect  of  which, 
unfortunately,  is  to  smear  the  finer  details  and  introduce  important 
inaccuracies  in  the  results.  These  effects  tend  to  be  reduced  as  the 
grid  resolution  increeses,  but  so  far  we  have  been  limited  by  our  Micro 
Vax  II  capabilities;  ebout  an  order  of  magnitude  increase  in  these, 
with  the  attendant  improvement  in  the  results,  is  expected  as  we  move 
into  our  new  multi-processor  computer. 

Figs.  2  and  3  show  results  obtained  for  one  test  case.  Gross 
instability  prevented  solution  for  20  KA,  apparently  due  to  very  low 
densities  near  the  anode.  Even  at  10  KA,  a  30k  reduction  in  anode 
density  can  be  seen  in  Fig.  2,  with  a  strong  density  gradient  on  the 
cathode  surface.  The  contours  in  Fig.  3  are  those  of  the  velocity 
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divided  by  the  local  aagnetoacoustic  speed,  which  theory  indicates 
should  take  the  place  of  the  ordinary  sound  speed  in  high-interaction 
flows. 

In  parallel  with  this  work,  we  are  currently  adapting  the  net  hod 
of  Flux-Corrected  Transport  (Ref.  1),  which  holds  greater  proaise  for' 
fully  three-diaensional  capabilities  when  coupled  with  the  new  vector 
aultiprocessor  aachine.  Preliminary  results  of  this  work  are  expected 
to  be  available  for  presentation  at  the  aeeting. 
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Year:  1976 

Laboratory:  Japan 

Publication:  AIAA  Paper  76-1002,  International  Electric  Propuls 
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Title:  MPD  Arcjet  Test  in  a  Large  Space  Chamber 


Year :  0 

Laboratory:  Japan 

Publication:  ISAS  Research  Note  33,  Institute  for  Space  and  Ast 
ronautical  Science,  Tokyo,  Japan 

Record  =  60 

Author:  Kuriki,  K: ,  'MaTsuo,  S. ,  Shoji,'  H. 

Title:  Metal  Plasma  Accelerator 


Year:  1979 

Laboratory:  Japan 

Publication:  AIAA  Paper  79-2092,  14th  International  Electric  Pr 

opulsion  Conference,  Princeton,  NJ,  Oct.  30-Nov.  1,  1979 
,  1979 

Record  =  61 

Author:  Kuriki,  K.,  Nakamaru,  K.,  Morimoto,  S. 

Title:  MPD  Thruster  Test  on  Engineering  Test  Satellite 


Year:  1979 

Laboratory:  Japan 

Publication:  AIAA  Paper  79-2071,  14th  International  Electric  Pr 

opulsion  Conference,  Princeton,  NJ,  Oct.  30-Nov.  1,  1979 
,  1979 

Record  =  62 

Author:  Kuriki,  K. 

Title:  The  MPD  Thruster  Test  on  the  Space  Shuttle 


Year :  1973 

Laboratory:  Japan 

Publication:  AIAA  Paper  79-661,  13th  International  Electric  Pro 

pulsion  Conference,  San  Diego,  CA,  April  25-27,  1973 
73 

Record  =  63 

Author:  Kuriki,  K.,  Onishi,  M.,  Morimoto,  S. 

Title:  Thrust  Measurement' of  K  III  MPD  Arcjet 


Year:  1981 

Laboratory:  Japan 

Publication:  AIAA  Paper  31-0*33,  15th  International  Electric  Pr 

opulsion  Conference,  April  21-23,  1931,  Las  Vegas,  NE 
NE 

Record  =  64 

Author:  Kunii,  Y.,  Shimizu,  Y.,  Kuriki,  K. 

Title:  Current  Distribution  in  a  Quasi-Steady  MPD  Arcjet 


Year:  1932 

Laboratory:  Japan 

Publication:  AIAA  Paper  32-1917,  16th  international  Electric  Pr 

opulsion  Conference,  New  Orleans,  LA,  Nov.  17-19,  1982 
1982 


Record  =65 

Author:  Kuriki,  K.,  Kunii,  Y.,  Shimizu,  Y. 

Title:  Idealized  Model  for  Plasma  Acceleration  in  an  MHD 
Channel 

Year:  1981 

Laboratory:  Japan 

Publication:  AIAA  Paper  81-0685,  15th  International  Electric  Pr 

opulsion  Conference,  Las  Vegas,  NE,  April  21-23,  1981 

981 

Record  =  66 

Author:  Kuriki,  K.,  Ijichi,  K.,  Harada,  H. ,  Okamura,  T. 

Title:  Mass  Reduction  of  Capacitor  Bank  for  MPD  Thruster 


j  Year:  1982 

Laboratory:  Japan 

Publication:  AIAA  Paper  82-1379,  16th  International  Electric  Pr 

opulsion  Conference,  New  Orleans,  LA,  Nov..  17-19,  1992 
1982 

Record  =  67 

■Author:  Ijichi,  K.,  Yoshida,  T.,  Kudo,  I.,  Kuriki,  K. 
j  Title:  Radiated  Emission  Noise  of  the  Plasma 


Year:  1982 

Laboratory:  Japan 

Publication:  AIAA  Paper  82-1883,  16th  International  Electric  Pr 

opulsion  Conference,  New  Orleans,  LA,  Nov.  17-19,  1982 
1982 

Record  =  68 

Author:  Kurtz,  H.L. ,  Auweter-Kurtz ,  M.,  Schrade,  H.O. 

Title:  Self-Field  MPD  Thruster  Design-  Experimental  and  T 
heoretical  Investigations 

Year :  1985 

Laboratory:  Germany 

Publication:  AIAA  Paper  85-2002,  13th  International  Electric  Pr 

opulsion  Conference,  Spet-  30-0ct.  2,  1985,  Alexandria 
dria,  VA 

Record  =  69 

Author:  Krulle,  G. 

Title:  Theoretical  Treatment  of  Current,  Mass,  Flow,  and 
Related  Distributions  in  MPD  Plumes 

Year:  1972 

Laboratory:  Germany 

Publication:  AIAA  Paper  72-501,  9th  Electric  Propulsion  Confere 
nee,  Bethesda,  MD,  April  17-19,  1972 

Record  =  70 

Author:  Hugel,  H. 

Title:  Self-Magnetic  Effect  in  Arcjet  Engines 


Year:  1967 

Laboratory:  Germany 

Publication:  AIAA  Journal  Technical  Notes,  Vol.  6,  No.  8,  Augus 
t,  1963 


Record  = 
Author : 
Title: 


71 

Hugel,  H.,  Kruelle,  G.,  Peters,  T. 

Investigations  on  Plasma  Thrusters  with  Thermal  an 
d  Self-Magnetic  Acceleration 


Year:  1967 

Laboratory:  Germany 

Publication:  AIAA  Journal,  Vol.  5,  No.  3,  March,  1967 


Record  =  72 

Author:  Schrade,  H.O.,  Auwetei — Kurtz,  M. ,  Kurtz,  H.L. 
Title:  Cathode  Erosion  Studies  on  MPD  Thrusters 


Year:  1985 

Laboratory:  Germany 

Publication:  AIAA  Paper  85-2019,  18th  International  Electric  Pr 

opulsion  Conference,  Sept.  30-0ct.  2,  1985,  Alexandria,  va 
dria,  VA 

Record  =  73 

Author:  Kurtz,  H.L.,  Auweter-Kurtz ,  M. ,  Schrade,  H.O. 

Title:  Experiments  with  a  Quasi-Steady  Coaxial  Ablation  M 
PD  Thruster  -  - 

Year :  1985 

Laboratory:  Germany 

Publication:  AIAA  Paper  85-2063,  19th  International  Electric  Pr 

opulsion  Conference,  Sept.  30-0cfc .  2,  1935,  Alexandria,  VA 
dria,  VA 

Record  =  74 

Author:  Schrade,  H.O.,  Auweter-Kur tz ,  M. ,  Kurtz,  H.L. 

Title:  Stability  Problems  in  MAgneto  Plasmadynamik  ARC  Th 
rusters 

Year:  1985 

Laboratory:  Germany 

Publication:  AIAA  Paper  35-1633,  13th  Fluid  Dynamics  and  Plasma 

dynamics  and  Lasers  Conference,  July  16-13,  1935,  Cincinnat 
Cincinnati,  OH 
Record  =  75 

Author:  Auweter-Kurtz ,  M. ,  Kurtz,  H.L.,  Schrade,  H.O. 

Title:  Optim.  of  Propulsion  Systems  for  Orbital  Transfer 

with  Separate  Powersupplies  Considering  variable  Thruster  Effs. 
hruster  Effs. 

Year:  1985 

Laboratory:  Germany 

Publication:  AIAA  Paper  35-1152,  2lst  Joint  Propulsion  Conferen 
ce,  July  3-10,  1935,  Monterey,  CA 

Record  =  76 

Author:  Schrade,  H.O. 

Title:  Comment  to  the  Stability  Theory  for  a  Current-Carr 
ying  Plasma  Channel 

Year:  0 

Laboratory:  Germany 

Publication:  No  publication  information  available. 


Record  =  77 

Author:  Krulle,  G.,  Zeyfang,  E. 

Title:  Preliminary  Conclusions  of  Continuous  Applied  Fiel 
d  Electromagnetic  Thruster  Research  at  DFVLR 


Year:  1975 

Laboratory:  Germany 

Publication:  AIAA  Paper  75-417,  llth  Electric  Propulsion  Confer 
ence,  New  Orleans,  LA,  March  19-21,  1975 


Record  =  78 

Author:  Buhler,  R.D. 

Title:  Plasma  Propulsion  for  Near-Earth  Missions  of  Large 
Space  Structures-  Status,  Problems,  and  Prospects 


Year:  1982 

Laboratory:  Germany 

Publication:  IRA-82  P7,  University  of  Stuttgart  Publication,  St 
uttgart,  FRG 


Record  =  79 

Author:  Buehler,  R.D.,  Lo,  R.E. 

Title:  Space  Propulsion  Systems.  Present  Performance  Lim 
its  and  Application  and  Development  Trends 


Year:  1981 

Laboratory:  Germany 

Publication:  Presented  at  German  Aerodyn.  &  Space  Flight  Congre 

ss,  Darmstad,  FRG,  Sept.  19-23,  1973.  Translation:  NASA  TM-761 
:  NASA  TM-761 63 
Record  =  30 

Author:  Auweter-Ming,  M. ,  Schrade,  H.O. 

Title:  Exploration  of  Arc-Spot  Motion  in  the  Presence  of 
Magnetic  Fields 


Year :  1980 

Laboratory:  Germany 

Publication:  Journal  of  Nuclear  Materials  93  &  *4,  (1930),  799- 
805 


Record  =  31 

Author:  Di  Capua,  M.S.,  Kurtz,  H.L. 

Title:  Electrolytic  Capacitor  Power  Source  Design  for  Qua 
si-Steady  Magnetoplasmadynamic  Arcs 


Year :  1973 

Laboratory:  Germany 

Publication:  AIAA  Paper  73-1107,  10th  Electric  Propulsion  Confe 
rence.  Lake  Tahoe,  NE,  Oct.  31-Nov.  2,  1973 


Record  =  32 

Author:  Guile,  A.E.,  Juttner,  B. 

Title:  Basic  Erosion  Processes  of  Oxidized  and  Clean  Meta 
1  Cathodes  by  Electric  Arcs 


Year:  1930 

Laboratory:.  Germany 

Publication:  IEEE  Transactions  on  Plasma  Science,  Vol.  PS-8,  No 
.  3,  September,  1980 


Record  =  33 

Author:  Baksht,  F.G.,  Moizhes,  B.Y.,  Rybakov,  A . B . 
Title:  Critical  Regime  of  a  PLasma  Accelerator 


Year:  1973 

Laboratory:  USSR 

Publication:  Sov.  Phys.  Tech.  Phys.,  Vol.  13,  No.  12,  June  1974 


Record  =  84 

Author:  Demidenko,  I. I.,  Krupnik,  L.I.,  Lomino,  N.S.,  Ovcharenko,  v. 
Title:  Operation  of  a  Quasisteady  Plasma  Accelerator  with 
a  Light-Gas  Anode  Sheath 

Year :  1979 

Laboratory:  USSR 

Publication:  Sov.  J.  Plasma  Physics,  Vol.  6,  No.  3,  May-June,  1 
980 

Record  =  05 

Author:  Glotova,  N.N.,  Kurdyashev,  K.P.,  Ostretsov,  I.N.,  Porotnikov 
Title:  RF  Waves  in  Electrodynamic  Plasma  Acceleration 


Year:  1975 

Laboratory:  USSR 

Publication:  Sov.  Phys.  Tech.  Phys.,  Vol.  21,  No.  3,  March,  197 
6 

Record  =  36 

Author:  Skoblik,  I.P.,  ZoLatotrubov ,  I.M.,  Novikov,  Y.M. 

Title;  Effect  of  Initial  Gas  Conditions  in  a  Coaxial  Acce 
lerator  on  the  Plasma  Parameters 

Year:  1972 

Laboratory:  USSR 

Publication:  Sov.  Phys.  Tech.  Phys.,  Vol.  13,  No.  2,  August,  19 
73 

Record  =  37 

Author:  Glotova,  Kazanskii,  Kirdyashev,  Ostretsov,  Porotnikov,  Utkin 
Title:  Anomalous  Operation  of  an  Elec trodynamic  Plasma  Ac 
celerator 

Year:  1977 

Laboratory:  USSR 

Publication:  Sov.  Phys.  Tech.  Phys.,  Vol.  23,  No.  7,  July,  1973 


Record  =  33 

Author:  Klimenko,  Kovalev,  Lykasov,  Lyapin,  Tkachev,  Chursin 
Title:  Influence  of  External  Magnetic  Field  on  the  Cathod 
e  Region  in  a  High-Current  Vacuum  Discharge 

Year:  1976 

Laboratory:  USSR 

Publication:  Sov.  Phys.  Tech.  Phys.,  Vol.  21,  No.  10,  October, 
1976 


Record  =  3'? 

Author:  Zhurin,  V.V.,  Porotnikov,  A. A.,  Porotnikov,  A. A.  Jr.,  Shadov 
Title:  MPD-Thruster ,  Opportunities,  and  Perspectives 


Year:  1981 

Laboratory:  USSR 

Publication:  AIAA  Paper  31-0689,  15th  International  Electric  Pr 

opulsion  Conference,  April  21-23,  1931,  Las-  Vegas,  NE 

NE 

Record  =  90 

Author:  Turchi,  P.J. 

Title:  The  Cathode  Region  of  a  Quasi-Steady  Magnetoplasma 
dynamic  Arcjet 

Year:  *  1970 

Laboratory:  Princeton 

Publication:  Ph.D.  Thesis  of  P.J.  Turchi,  Department  of  Aerospa 
ce  and  Mechanical  Sciences,  Princeton  University 

Record  =  91 

Author:  Tang,  M.,  Jacobson,  D. 

Title:  Electron  Emission  from  Cesiated  IX  and  2X  Thoriate 
d  Tungsten  Alloys 


Year:  1994 

Laboratory:  Arizona 

Publication:  AIAA  Paper  84-1759,  19th  Thermophysics  Conference, 
June  25-23,  1934,  Snowmass,  CO 

Record  =  92 

Author:  Kimblin,  C.W. 

Title:  Cathode  Spot  Erosion  and  Ionization  Phenomena  in  t 

he  Transition  from  Vacuum  to  Atmospheric  Pressure  Arcs 
Arcs 

Year:  1974 

Laboratory:  Westinghouse 

Publication:  Journal  of  Applied  Physics,  Vol .  45,  No.  12,  Decern 
ber,  1974 

Record  =  93 

Author:  Lutz,  M . A . 

Title:  Radiation  and  its  Effect  on  the  Non-Equilibrium  Pr 
operties  of  a  Seeded  Plasma 

Year :  1965 

Laboratory:  MIT 

Publication:  Ph.D.  Thesis  of  m.a.  Lutz,  Department  of  Aeronauti 
cs  and  Astronautics ,  M.I.T. 

Record  =  94 

Author.-  Burton,  R.L.,  Goldstein,  S.A.,  Hilko,  B.K.,  Tidman,  D.A.,  wi 
Title:  Investigation  of  a  Pulsed  Elec trothermal  Thruster 


Year:  1983 

Laboratory:  GT  Devices 

Publication:  NASA  CR-16826*,  GTD  83-10,  NASA  Lewis  Research  Cen 
ter  Contract.  NAS  3-2377-3 


■>  ->  •> 


Record  =  95 

Author:  Michels,  C.J.,  York,  T.M. 

Title:  Pressure  Measurements  in  the  Exhaust  of  a  Pulsed  M 
egawatt  MPD  Arc  Thruster 

Year:  1971 

Laboratory:  NASA  Lewis 

Publication:  AIAA  Paper  71-196,  9th  Aerospace  Sciences  Meeting, 
New  York,  NY,  January  25-27,  1971 

Record  =  96 

Author:  Michels,  C.J. 

Title:  Dynamic  Voltage-Current  Charac teristics  of  a  Megaw 
att  MPD-ARC  Thruster 

Year:  1970 

Laboratory:  NASA  Lewis 

Publication:  AIAA  Paper  70-164,  3th  Aerospace  Siences  Meeting, 
New  York,  NY,  January  19-21,  1970 

Record  =  97 

Author:  John,  R.R.,  Bennett,  S.,  Connors,  J.F. 

Title:  Arcjet  Engine  Performance:  Experiment  and  Theory 


Year:  1963 

Laboratory:  Avco 

Publication:  AIAA  Journal,  November,  1963  , 


Record  =  90 

Author:  Choueiri,  E.Y. 

Title:  A  Scaling  Strategy  for  the  Preliminary  Design  of  M 
agnetoplasmadynamic  Thrusters 

Year:  1982 

Laboratory.-  Syracuse  U. 

Publication:  Master’s  Thesis,  Syracuse  University,  December,  19 
83 

Record  =  99 

Author:  Rudolf,  L.K.,  Hamlyn,  K.M.,  Ogg,  G.M.,  Davis,  H.P.,  Stump,  W 
Title:  MPD  Thruster  Definition  Study 


Year:  1984 

Laboratory:  AFRPL 

Publication:  Air  Force  Rocket  Propulsion  Laboratory  Publication 
AFRPL  TR-34-046 

Record  =  100 

Author:  Poon,  T.W. 

Title:  Local  Stability  Anaysis  of  the  Magnetoplasmadynami 
c  (MPD)  Thruster 

Year:  1985 

Laboratory;  MIT 

Publication:  Master’s  Thesis,  Department  of  Aeronautics  and  Ast 
ronautics,  M.I.T.,  September,  1935 


Record  =  101 

Author:  Malliaris,  A.C.,  John,  R.R.t  Garrison,  R.L.,  Libby,  D.R. 
Title:  Performance  of  Quasi-Steady  MPD  Thrusters  at  High 
Powers 

/ 

Year :  1971 

Laboratory:  Avco 

Publication:  AIAA  Journal,  Vol.  10,  No.  2,  February,  1972.  Als 
O  N71 -33543. 

Record  =  102 

Author:  Patrick,  R.M.,  Schneiderman,  A.M. 

Title:  Performance  Character istics  of  a  Magnetic  Annular 
Arc 

Year:  1966 

Laboratory:  Avco 

Publication:  AIAA  Journal,  Vol.  4,  No.  2,  February,  1966. 


Record  =  103 

Author:  Schneiderman,  A.M.,  Patrick,  R.M. 

Title:  Optimization  of  the  thermal  Efficiency  of  the  Magn 
etic  Annular  Arc 

Year:  1966 

Laboratory:  Avco 

Publication:  AIAA  Journal,  vol.  4,  No.  10,  October,  1966. 


Record  =  104 

Author:  Chubb,  D.L. 

Title:  Fully  Ionized  Guasi-One  Dimensional  Magnetic  Nozzl 
e  Flow 

Year:  1971 

Laboratory:  NASA  Lewis 

Publication:  AIAA  Paper  71-141,  9th  Aerospace  Sciences  Meeting, 
New  York,  NY,  January  25-27,  1971 

Record  =  105 

Author:  John,  R.R.,  Bennett,  S.,  Connors,  J.F. 

Title:  Experimental  Performance  of  a  High  Specific  Impuls 
e  Arc  Jet  Engine 

Year:  1965 

Laboratory:  Avco 

Publication:  Astronautica  ACTA,  Vol.  11,  No.  2,  1965 


Record  =  106 

Author:  Koester,  J.K.,  Dick,  R.S.,  Britt.,  E.J. 

Title:  Operating  Character istics  of  Space  Nuclear  Reactor 
s  Coupled  to  Electric  Propulsion  Spacecraft 

Year:  1985 

Laboratory:  Space  Power 

Publication:  Presented  at  IEEE  Conference  on  Plasma  Science  and 
Workshop'  on  MPD  Thrusters,  June  3-4,  1985 


Author:  Wetch,  J.R.,  Koester,  J.K.,  Britt,  E.J. 

Title:  The  Safety  Impact  of  NEP  From  LEO  to  Operational  0 
rbits 

Year :  1985 

Laboratory:  Space  Power 

Publication:  Presented  at  Workshop  on  Systems  and  Fundamental  R 

esearch  of  MPD  Thrusters,  June,  1985,  Pittsburgh,  PA 

PA 

Record  =  108 

Author:  John,  R.R.,  Bennett,  S. 

Title:  Recent  Advances  in  Electrothermal  and  Hybrid  Elect 
rothermal-Electromagnetic  Propulsion 

Year:  0 

Laboratory:  Avco 

Publication:  Publication  Information  Unavailable 


Record  =  109 

Author:  Lawless,  J.L. 

Title:  Efficiency  and  Onset  in  MPD  Thrusters 


Year :  1983 

Laboratory:  CMU 

Publication:  Department  of  Mechanical  Engineering  Publication, 
Carnegie-Mellon  University,  Pittsburgh,  PA 

Record  =  110 

Author:  Turchi,  P.J. 

Title:  Application  of  Hall  Effect  Inductive  Switching  to 
Quasi-Steady  MPD  Arcjets 

Year:  1979 

Laboratory:  NRL 

Publication:  AIAA  Paper  79-2085,  14th  International  Electric  Pr 

opulsion  Conference,  Princeton,  NJ,  Oct.  30-Nov.  1,  1979 
,  1979 

Record  =  111 

Author:  Martinez-Sanchez ,  M. ,  Heimerdinger ,  D. 

Title:  Two-Dimensional  Analysis  of  an  MPD  Arcjet 

Year:  1985 

Laboratory:  MIT 

Publication:  AIAA  Paper  35-2040,  13th  International  Electric  Pr 

opulsion  Conference,  Sept.  30-0ct.  2,  1985,  Alexandria,  VA 
dria,  VA 

Record  =  112 

Author.-  Orlandi,  A.,  Baiocchi,  G.  ,  Carabella,  G.,  Andrenucci,  M. 

Title:  Solid  Propellant  MPD  Thruster  Development 


Year:  1985 

Laboratory:  Italy 

Publication:  AIAA  Paper  85-2066,  18th  International  Electric  Pr 

opulsion  Conference,  Sept.  30-0ct.  2,  1935,  Alexandria,  VA 
dria,  VA 


Record  =  113 

‘Author:  Andrei,  G.,  Andrenucci,  M.,  Choaeiri,  E.Y. 
Title:  Theoretical  Aspects  of  MPD  Thruster  Scaling 


Year :  1984 

Laboratory:  Italy 

Publication:  IEPC  84-10,  17th  International  Electric  Propulsion 
Conference,  Tokyo,  1984 

Record  =  114 

Author:  Lawless,  J.L.,  Subramanian,  V.V. 

Title:  A  Theory  of  Onset  in  MPD  Thrusters 


Year:  1983 

Laboratory:  CMU 

Publication:  Publication  of  Carnegie-Mellon  University,  Pittsbu 
rgh ,  PA 

Record  =  115 

Author:  Michels,  C.J.,  Rose,  J.R.,  Sigman,  D.R. 

Title:  Electron  Number  Density  and  Temperature  in  MPD-Arc 
Thruster  Exhausts 

Year :  1972 

Laboratory:  NASA  Lewis 

Publication:  AIAA  Paper  72-209,  10th  Aerospace  Sciences  Meeting 
,  San  Diego,  CA,  January  17-19,  1972 

Record  =  116 

Author:  Rudolph,  L.K.,  King,  D.Q. 

Title:  Electromagnetic  Thrusters  for  Spacecraft  Prime  Pro 
pulsion 

Year:  1984 

Laboratory:  JPL 

Publication:  AIAA  Paper  84-1446,  20t.h  Joint  Propulsion  Conferen 
ce,  June  11-13,  1934,  Cincinnati,  OH 

Record  =  117 

Author:  Heimerdinger ,  D.J. 

Title:  An  Approximate  Two-Dimensional  Analysis  of  an  MPD 
Thruster 

Year:  1984 

Laboratory:  MIT 

Publication:  Master’s  Thesis,  Department  of  Aeronautics  and  Ast 
ronautics,  M.I.T.,  June,  1984 

Record  =  118 

Author:  Lawless,  J.L.,  Subramanian,  V.V. 

Title:  A  Theory  of  Onset  in  Magnetoplasmadynamic  Thruster 
s 


Year :  1935 

Laboratory:  CMU 

Publication:  AIAA  Paper  85-2039,  13th  International  Electric  Pr 

opulsion  Conference,  September  30-0ctober  2,  1985,  Alexandria, 
Alexandria,  va 


.Record  =  119 

Author:  Rock,  B.A.,  Mantenieks,  M.a.  ,  Parsons,  M.L. 

Title:  Rapid  Evaluation  of  Ion  Thruster  Lifetime  Using  Op 
tical  Emission  Spectroscopy 

Year :  1985 

Laboratory:  Arizona,  NASA  Lewis 

Publication:  AIAA  Paper  95-2011,  13th  International  Electric  Pr 

opulsion  Conference,  September  30-0ctober  2,  1985,  Alexandr 
Alexandria,  VA 

Record  =  120 

Author:  Michels,  C.J.,  Sigman,  D.R. 

Title:  Exhaust  Character istics  of  a  Megawatt  Nitrogen  MPD 
-ARC  Thruster 

Year :  1971 

Laboratory:  NASA  Lewis 

Publication:  AIAA  Journal,  Vol.  9,  No.  6,  June  1971 


Record  =  121 

Author:  Cox,  J.E. 

Title:  Electromagnetic  Propulsion  Without  Ionization 


Year:  1931 

Laboratory : 

Publication:  Journal  of  Spacecraft  and  Rockets,  Vol.  18,  No.  5,  . 
September-October  1931. 

Record  =  122 

Author:  Wolff,  M. ,  Kelly,  A. J. ,  Jahn ,  R.G. 

Title:  A  High  Performance  Magnetoplasmadynamic  Thruster 


Year:  1934 

Laboratory:  Princeton 
Publication:  IEPC  34-32 


Record  =  123 

Author:  Mehta,  R.C. 

Title:  Heat  Transfer  Analysis  of  a  Hollow  Cathode  of  an  M 
PD  Arc  Thruster 

Year:  1984 

Laboratory:  India 
Publication:  IEPC  34-12 


Record  =  124 

Author:  Cassady,  R.J.,  Stoyanof,  D.S.,  Kilgore,  O.W. 

Title:  Characterization  of  Performance  variation  With  Geo 
metry  in  an  MPD  Thruster 

Year :  1935 

Laboratory:  AFRPL 

Publication:  AIAA  Paper  85-2021,  18th  International  Electric  Pr 

opulsion,  September  30-0ctaber  2,  1935,  Alexandria,  VA 
,  VA 


'  Record  =  125 

Author:  Brushlinsky,  K.V.,  Kozlov,  A.N. 

Title:  Numerical  Simulations  of  Self-Ionizing  Gas  Flows 


Year :  1983 

Laboratory:  USSR 
Publication:  AIAA  Paper  33-1916 


Record  =  126 

-Acrttfar  r*Tarrara  ,■  tt.  T~*Tmu  r'a'rrr -  ' 

Title:  Theoretical  Analysis  of  Current  Distribution  and  P 
lasma  Acceleration  in  MPD  Arcjet 

l 

nYear :  1934 

Laboratory:  Japan  (Tokyo) 

Publication:  IEPC  34-07 


Record  =  127 

Author:  Turchi,  P.J.,  Boyer,  C.N.,  Davis,  J.F. 
Title:  Multi-Stage  Plasma  Propulsion 


Year:  1984 

Laboratory:  RDA 

Publication:  IEPC  34-51  JSASS/AIAA/DGLR  17th  International  Elec 
trie  Propulsion  Conference,  Tokyo,  Japan,  1984 

Record  =  123 

Author:  Maisenhalder ,  F. ,  Mayerhofer,  W. 

Title:  Jet-Diagnostics  of  a  Self-Field  Accelerator  with  L 
angmuir  Probes 

Year:  1974 

Laboratory:  Germany  (Stuttgart) 

Publication:  AIAA  Journal  Vol.  12,  No.  9.  September  1974 


End  of  locate  scope 


